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Abstract 

First, we show that Sturm algorithm and Sylvester algorithm, which compute the number 
of real roots of a given univariate polynomial, lead to two dual tridiagonal determinantal 
representations of the polynomial. Next, we show that the number of real roots of a 
polynomial given by a tridiagonal determinantal representation is greater than the signature 
of this representation. 

Introduction 

There are several methods to count the number of real roots of an univariate polynomial 
p(x) £ M.[x] of degree n (for details we refer to [BPR]). Among them, the Sturm algorithm 
says that the number of real roots of p(x) is equal to the number of Permanence minus the 
number of variations of signs which appears in the leading coefficients of the signed remainders 
sequence of p(x) and p'(x). 

Another method is the Sylvester algorithm which says that the number of real roots of p(x) 
is equal to the signature of the symmetric matrix whose (i,j)-th entry is the i + j-th Newton 
sums of the roots of the polynomial p(x). 

One purpose of the paper is to point out, at least in the generic situation, that these two 
classical algorithms can be viewed as dual. 

In section 1, we introduce signed remainders sequences of two given monic polynomials p(x) 
and q(x) of respective degrees n and n — 1. With some conventions of signs and others, we 
give a presentation of this sequence through a tridiagonal matrix Td(p, q). Next, we give a 
decomposition of this tridiagonal matrix as Td(p, q) = LC^L" 1 where L is lower triangular and 
CT is the transposed of the companion matrix associated to p(x). 



In section 2, we introduce the duality between the Sturm and Sylvester algorithm, first when 
the polynomial p(x) has only single and real roots, and then in Theorem 2.5 we generalize it to 
the generic case. 

More precisely, on one hand we have 

f p(x) = det(xld n - Td(p, q)) 

\ q(x) = det(xld n _i - Td(p,g) n _i) 

with the conventions that Id„ (or Id in short) denotes the identity matrix of M nxn and 6 R kxk 
(respectively A^ G R kxk ) denotes the k-th principal submatrix (respectively the k-th antiprincipal 
submatrix) of A which corresponds to extracting the first k (respectively the last k) rows and 
columns in the matrix A £ ]R nxn . 

On the other hand, we consider a natural Hankel (hence symmetric) matrix H(q/p) 6 M nxn 
associated to p(x) and q{x). Generically it admits an LU decomposition of the form H(q/p) = 
KJK T where J is a signature matrix (a diagonal matrix with coefficients ±1 onto the diagonal) 
and K is lower triangular. Then, we introduce the tridiagonal matrix Td = K _1 Cp K , which is 
such that p{x) = det(xld n — Td). 

If we consider that the matrices Td(p, q) and Td represent linear mappings in some basis, 
then the duality Theorem 2.5 means that one matrix can be deduced from the other simply by 
reversing the ordering of the basis. 

We shall mention that, in the case when all the roots of p{x) are real, the existence 
of a tridiagonal symmetric matrix Td given by the signed remainders sequence of p(x) and 
q{x) together with the identity p{x) = det(xld n — Td) corresponds to the Routh-Lanczos 
algorithm which answers a structured Jacobi inverse problem. Namely, the question to find a 
real symmetric tridiagonal matrix A with a given characteristic polynomial p{x) and such that 
the characteristic polynomial of its principal minor A n -i, of size n — 1, is proportional to p'{x). 
We refer to [EP] for a survey on the subject. 

In section 3, we focus on the relation of the question of real roots counting and the question of 
determinantal representation. We say that p(x) = det( J— xA) is a determinantal representation 
of the polynomial p(x) if J £ R nxn is a signature matrix and A 6 M. nxn is a symmetric matrix. 

Remark that we may transform the identity p(x) = det(J — xA) into p*{x) = det(xJ — A) 
where p*(x) is the reciprocal polynomial of p(x). If we write 

p*(x) = det(J) x det(zld- AJ), 

it shows a connexion with the results of section 2 when the matrix AJ is tridiagonal. More 
precisely, we establish that such a determinantal representation is always possible : we may 
even find a family of representations for a given polynomial p(x). We show also that, given such 
a determinantal representation for a polynomial p(x), its number of real roots is at least equal 
to the signature of the signature matrix J. 

Finally, in section 4 we end with some worked examples. 
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1 Tridiagonal representation of signed remainders sequences 
1.1 Definitions 

Let a = (cci, . . . , a n ), (3 = (Pi, . . . , (3 n -i) an d 7 = (71 j ■ • • , 7n-i) be three sequences of real 
numbers. We set the tridiagonal matrix Td(a, j3, 7) to be : 

/ a n 7„_i ... \ 

Pn-i 0i n -\ 7„_ 2 : 

Td(a,/3,7)= /3 n _ 2 "•■ "■■ 

: '•• '•• 71 

\ ... ft ai y 

Let p(x) and g(x) be two monic polynomials of respective degrees n and n — 1. We set 
SRemS(p, g) = (p k (x)) k to be the signed remainders sequence of and g(x) defined in the 
following way : 



(1) 



where 



Po(x) = p(x) 
pi(x) = q(x) 
Pk{x) = q k+ i(x)p k+1 (x) 



e k +iPl +1 Pk+2{x) 



(2) 



f Pk{x),q k+1 (x) £ I 
e k+1 G {-!,+!}, 



x 



(3 k +i is a positive real number, 
p k+ 2(x) is monic and degp fc+ 2 < degp k+ i. 

This is a finite sequence which stops at the step just before we reach the zero polynomial as 
remainder. With these conventions, the signed remainders sequence SRemS(j>, q) that we obtain 
is also called the Sturm-Habicht sequence of p(x) and q(x). 

Let us assume that there is no degree breakdown in SRemS(p, q). Namely : 
(3) (Vfc€{0,...,n}) (deg Pk = n-k) 

Then, q k+ i(x) is a degree one polynomial which we write q k+ \(x) = (x — a k+ i) with a k+ \ G M. 
Another consequence is that gcd(p, q) = 1. 

Let 7fc+i = e k+ \(5 k+ \ and consider the following tridiagonal matrix : 

Td(p, g ) = Td(a,/?,7) 

We may read on this matrix all the informations about the signed remainders sequence 
SRemS(p, q). 

For a given tridiagonal matrix Td = Td(a,/3, 7) G M nxn , we define the /irsi principal 
lower diagonal (respectively the first principal upper diagonal) of Td to be the sequence 
P = (Pi, . . . , (3 n -i) (respectively 7 = (71, . . . , 7 n _i)). We will say that these first principal 
diagonals are non-singular if all the coefficients Pi (respectively 7,) are different from zero. 

Note that the no degree breakdown assumption (3) implies that the principal diagonals of 
Td(p, q) are non-singular. 
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Proposition 1.1 (i) To any tridiagonal matrix Td = Td(a,/3, 7) with non- singular principal 
diagonals, we may canonically associate a (unique) couple of monic polynomials p(x) and 
q{x) of respective degrees n and n — 1 such that the sequence SRemS(p, q) has no degree 
breakdown and the characteristic polynomial o/Tdfc is equal to p n -k( x ) ■' 

det(xld fc - Td fc ) = p n -k(x). 

(ii) To any couple of monic polynomials p(x) and q(x) of respective degrees n and n — 1 such 
that SRemS(p, q) has no degree breakdown, we may associate a unique tridiagonal matrix 
with non-singular principal diagonals Td(p, q) = Td(a,{3, 7) satisfying for all k, (5^ > 
and 7fc = ekfik where e& = ±1. 

(Hi) When we have (i) and {ii), the matrix Td(p,q)) x P is tridiagonal and symmetric, where 
we have set 

I e n _i x...xei \ 
P = e 2 x ei 

V 1 / 

(iv) When we have (i) and (ii), the sequence of signs in the leading coefficients of the signed 
remainders sequence SRemS(p, q) is : 

(1, 1; ei, e 2, ei x £3, e 2 x £4, ei x e 3 x e 5 , . . . , e n _ 1 mo( j L 2 x . . . x e n _3 x e n _i) 

Proof : Concerning (i), the polynomials p(x) and q(x) are taken to be p(x) = det(xld n — Td) 
and q(x) = det(xld n _i — Td n _i). Then, we set for all k, 

(5 n _ fc (x) = det(xld fc - Td k ) 

(where Td^ is the k-th principal submatrix of Td) and we develop the determinant 

5 (x) = det(xld„ - Td) 

with respect to the last row. We get 

S (x) = (x - ai)(5i(x) - (/3i7i)<5 2 (x) 

Repeating the process, we obtain the same recurrence relation as the one defining the sequence 
(Pk( x ))k m (!)• Since Sq(x) = po(x) and <5i(x) = pi(x), we get the wanted identity. 

Point (ii) follows straightforward from the beginning of the section, whereas points (Hi) and 
(iv) follows from elementary computation. O 

We may note that to the tridiagonal matrix Td(p,q), we may associate also another natural 
polynomial remainder sequence : SRemS(j>, q) = SRemS(p, q) where 

p(x) = det(xld n - Td) 

and 

q(x) = det(xld n _i - Td n _i), 
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with the convention that Td& is the k-th. antiprincipal submatrix of Td. The signed remainders 
sequence SRemS(p, q) will be considered as the dual signed remainders sequence of SRemS(p, q). 
This only means that we may read on a tridiagonal matrix from the top left rather than from 
the bottom right ! 



For cosmetic reasons we will write Td(p, q) in place of Td(p, q). We obviously have : 



(4) Td(p, g)=Adx Td(p, q) x Ad 

where Ad„ G R nxn (Ad in short) stand for the anti-identity matrix of size n : 

/ ... 1 \ 

• •' •■' 



Ad„ = 



••' •■' : 
\l ... Oj 



1.2 Companion matrix 

We denote by A the transposed of the matrix A £ 
of the polynomial p(x) = x n + a n -\x 



C p — 



( 
1 


v° 

We recall a well-know identity (see for 



nxn and we define the companion matrix 
+ ao to be 

... -a \ 
: -a x 



'■■ 
1 



-a n -2 
-a n -i ) 



instance [EP]) : 



Proposition 1.2 Let p(x) and q(x) be two monic polynomials of respective degrees n and n—1 
such that SRemS(p, q) has no degree breakdown. 

Then, there is a lower triangular matrix L such that 

(5) Td(p,g) = LC p r L- 1 

Proof : With the notation of Subsection 1.1, let V{x) = (71 . . . 7„_ip n (x), . . . , 71^2(^)5 Pi (%))• 
A direct computation gives 

V{x) (Td(p, q)) T = xV(x) + (0, . . . , 0, -p(x)) 

Let U be the upper triangular matrix whose columns are the coefficients of the polynomials of 
V(x) in the canonical basis C(x) = (1, x, . . . , x n ~ l ). In other words : 



Besides, we have 



C(x)U = V(x) 



C{x)C p = xC(x) + (0, . . . , 0, -p(x)) 
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Thus 



C(x)C p U = xC(x)U + (0, . . . ,0, — p(x)) since pi(x) is monic 
= V(x)(Td(p,q)f 



= C(x)U(Td(p,q)) T 
We deduce the identity 

V( Xl , x n )C p U = V( Xl , x n )U (Td(p, q)f 



for any Vandermonde matrix V(xi, . . . , x n ) whose lines are (1, Xi, . . . , x™ ) for i = 1 . . . n. If we 
choose the n reals x±, . . . , x n to be distinct, then V(x±, . . . , x n ) becomes invertible and we get : 

Td(p,q)=LC^L- 1 



where L is the lower triangular matrix defined by L = U T . 



□ 



The following result says that the decomposition generically exists for any tridiagonal matrix, 
and also it is unique : 

Proposition 1.3 Any tridiagonal matrix Td with non-singular principal diagonals can be 
written Td = LCpL^ 1 where p{x) = det(xld — Td) and L is a lower triangular matrix. 
Moreover the matrix L is unique up to a multiplication by a real number. 

Proof : The existence is given by Proposition 1.1 and Proposition 1.2. 

We come now to the unicity. Assume that LiCpL^ 1 = L2CJL2 1 where L\ and L2 are lower 
triangular. Then, L = L^ l L\ is a lower triangular matrix which commute with Cj. 

If L = {ti,j)i<i,j< n , then 



LCl 



( t M 



in-1,1 
? 





tn— l,n— 1 
? 



and 



C T V L 



( t 2 ,i t 2 , 2 

*3,1 *3,2 *3,3 



\ 



tn,l 
? 





tn,n— 1 in,n 
? 



Thus ti,i = t 2 ,2 = . . . = t n:n and t 2 ,i = 1^3,2 = . . . = t n , n -i = and t 3j i = t 4j2 = . . . = 
tn n-2 = 0, and so on until t n 1 = 0. We deduce that L = Aid and we are done. O 
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1.3 Sturm algorithm 

As a particularly important case of signed remainders sequences, we shall mention the Sturm 
sequence which is SRemS(p, q) where q is taken to be the derivative of the polynomial p(x) up 
to normalization, i.e. q = p'/deg (p). 

For a given finite sequence v = {y\, . . . , UjS) of elements in { — 1, +1}, we recall the Permanence 
minus Variations number : 

fc-i 

PmV(i/i, . . . , v k ) = ^ Vin+i- 
i=i 

Here the sequence v will be for the sequence of signs of leading coefficients in SRemS(p, q). Then, 
the Sturm Theorem [BPR, Theorem 2.50] says that the number PmV(f) is exactly the number 
of real roots of p(x). 

If we assume that the polynomial p(x) has n distinct real roots, then the Sturm sequence 
has no degree breakdown and for all k we have = 1. Hence we get a symmetric tridiagonal 
matrix Td(p, q) which has the decomposition Td(p, q) = LCpL" 1 where L is the lower triangular 
matrix defined as in subsection 1.2. In particular, the last row of L gives the list of coefficients 
of the polynomial q(x) in the canonical basis. 

2 Duality between Sturm and Sylvester algorithms 
2.1 Sylvester algorithm 

Let us introduce the symmetric matrix Newt p (n) = {n>i,j)o<i,j<n-i define as 

n itj = Trace (C l p +j ) = N i+j 

which is nothing but the i + j-th Newton sum of the polynomial p(x). To be more explicit, if 
ai, . . . ,a n denote all the complex roots of the polynomial p(x), then the k-th Newton sum is 
the real number 7V fe = a\ + . . . + a^. 

Recall that the signature sign(^4) of a real symmetric matrix A G R nxn , is defined to be the 
number p — q, where p is the number of positive eigenvalues of A (counted with multiplicity) and 
q the number of negative eigenvalues of A (counted with multiplicity). The Sylvester Theorem 
(which has been generalized later by Hermite : [BPR, Theorem 4.57]) says that the matrix 
Newtp(n) is invertible if and only if p(x) has only single roots, and also that sgn(Newt p (n)) is 
exactly the number of distinct real roots of p{x). 

In particular, if the polynomial p(x) has n distinct real roots, then the matrix Newt p (n) is 
positive definite. Thus, by the Choleski decomposition algorithm, we can find a lower triangular 
matrix K such that Newt p (n) = KK T . Le us show how to exploit this decomposition. 

First, we write 

p{x) = det(xld - Cp) 

Then, we introduce a useful identity (which will be discussed in more details in the forthcoming 
section) : 

Newtp (n)C p = CjNewt p (n), 

So, we get : 

p(x) = det(xld - K^C^K) 
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Note that the matrix K Cj } K is tridiagonal. Our purpose in the following is to establish a 
connexion with the identity 

p{x) = det(xld - LC T p L- v ) 

obtained in Proposition 1.3. 

More generally, we will point out a connexion between tridiagonal representations associated 
to signed remainders sequences on one hand, and tridiagonal representations derived from 
decompositions of some Hankel matrices on the other hand. 



2.2 Hankel matrices and Intertwinning relation 

Roughly speaking, the idea of previous section is to start with the canonical companion identity 

p(x) = det(xld - Cp) 

and then to use a symmetric invertible matrix H satisfying the so-called intertwinning relation 



(6) 



HC P — CpH 



Since H is supposed to be symmetric invertible, Equation (6) only says that the matrix HC P 
is symmetric. It is a classical and elementary result that a matrix H satisfying equation (6) is 
necessarily an Hankel matrix. 



Definition 2.1 We say that the matrix H = (hij)o<ij< n -i G 



is an Hankel matrix if 



h 



1,3 



htfji whenever i + j = %' + f . Then, it makes sense to introduce the real numbers 
di+j = hi j which allow to write in short H = (ai+j)o<i > j< n -i- 

Let s = (sfc) be a sequence of real numbers. We denote by H n (s) or by H(sq, . . . , S2 n -2) the 
following Hankel matrix o/K riX ™ : 



H n {s) — (si+j)o<i,j<n-l 



( -s 

Si 



Si 



\ S n S n -)-l 



Sn \ 
Sn+1 

S2n-2 / 



We get from [BPR, Theorem 9.17] : 

Proposition 2.2 Let p(x) = x n + a^-ix™ -1 + . . . + oq and s = (s^) be a sequence of real 
numbers. The following assertions are equivalent 

(i) (Vfc > n) (s k = -a n _is fe _i - ... - a s k _ n ) 

(ii) There is a polynomial q(x) of degree degg < degp such that 

(Hi) There is an integer r < n such that det(H r (s)) / 0, and for all k > r, det(H k (s)) = 0. 
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Whenever these conditions are fulfilled, we denote by H n (q/p) the Hankel matrix H n (s). 

Back to the intertwinning relation (6) : it is immediate that an Hankel matrix H is a 
solution if and only if the (finite) sequence (so, . . . , §2^-2) satisfies the linear recurrence relation 
of Proposition 2.2(«), for k = n, . . . , 2n — 2. 

For further details and developments about the intertwinning relation, we refer to [HV]. 

The vector subspace of Hankel matrices in K nxn satisfying relation (6) has dimension n, and 
contains a remarkable element : the Hankel matrix Newt p (n) that was considered in subsection 
2.1 about Sylvester algorithm. Indeed, it is a well-known and elementary fact that the N^'s are 
real numbers which verify the Newton identities : 

(Vfc > n) (N k + a n _iiV fc _i + . . . + a N k _ n = 0) 



2.3 Barnett formula 

First, recall that if p(x) = x n + a Ti _ix n ~ 1 + . . . + ao and q(x) is a (non-necessarily monic) 
polynomial in R[x] whose degree is equal to n — 1, the Bezoutian of p(x) and q(x) is defined as 
the two-variables polynomial : 

Bez(p >g ) = g(y)p(a) - g(g)p(y) €R[x,y] 
x-y 

Let B(z) be any basis of the n-dimentional vector space WL[z]/p(z) over R. We denote by 
Bezg(p,q) the symmetric matrix of the coefficients of Bez(p,q) with respect to the basis B{x) 
and B(y). 

Among all the basis of R[z]/p(z) that will be interesting for the following, let us mention 
the canonical basis C = (l,z,..., z n ~ l ) and also the (degree decreasing) Horner basis TL{z) = 
(ho, . . . , h n _i) associated to the polynomial p(z) and which is defined by : 

' h (z) = z n ~ l + a n ^z n - 2 + . . . + 01 

hi(z) = z n - 1 - 1 + a n -xz n ^ % + ■ ■ ■ + a i+ i = zh i+1 (z) + a i+1 

h n -2(z) = z + a n -i 
k h n -i(z) = 1 



We recall from [BPR, Proposition 9.20] : 

Proposition 2.3 Let p(x) and q(x) be two polynomials such that degg < degp = n. Let s be 
the sequence of real numbers defined by 

p(x) 

Then, Bez n (p,q) = H n (s) = H n (q/p). 
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We come to a central proposition which is a consequence of the Barnett formula [Ba]. 

Proposition 2.4 Let p(x) and q(x) be two polynomials such that degg < degp = n and let Pqh 

be the change of basis matrix from the canonical basis C to the Horner basis H . We have 

q{C p ) = Pl n X H n (q/p) 

Proof : The Barnett formula has been established in [Ba] using direct matrix computations. 
For the convenience of the reader, we give here another proof (which may be found at various 
places in the literature). 
The obvious identity 

q(y)(p(x) - p(y)) = q(y)p(x) - p(y)q(x) + p(y)(q(x) - q(y)) 

implies, by definition of the Bezoutian B(p, q), that : 

/ \P( x ) ~ p(y) ~r~k / \ . > 

q(y) = Bez(p,q) mod p{y) 

x y 

Noticing that S^Kl = ^ hj ( y ) x J , we get 

n— 1 

q(y) h j(y) x3 = c (y)^ ez c(p, q)C(x) T mod p(y) 

3=0 

In other words, if we denote by M the matrix whose columns are the coefficients of q(y)hj(y) 
in the basis C(y), we get the identity 

C(y)MC{x) T = C(y)Bez c (p, q)C(x) T mod p{y) 

Since C p is the matrix of the multiplication by y mod p(y) with respect to the canonical 
basis C(y), we have also the identity 

M = q(C p )P Cn 

where the change of basis matrix Pen IS m f ac t the following Hankel matrix 

Pen = H(a ,ai,..., a n _i, 1,0,..., 0) G M nxn 
with the usual notation p{x) = x n + a n -\x n " 1 + . . . + ao- Hence, we get the Barnett Formula 

Bez c (p,q) = q(C p )P CH 
Finally, by Proposition 2.3, we derive the wanted relation : 

Bez c (p,q) = P? n Bez H (p,q)Pc H = Pl H H n {q/p)P CH 
Which concludes the proof. O 

To end the section, we show how Sturm and Sylvester algorithms can be considered as 
dual, in the case where all the roots of p{x) are real and simple, say x\ < . . . < x n . Then, 
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q(x) = p'(x)/n has also n — 1 simple real roots y± < . . . < y n ~\ which are interlacing those of 
p(x). Namely 

x 1 < yi < X 2 < V2 < ■ ■ ■ < Vn-l < x n 

We may repeat the argument to see that this interlacing property of real roots remains for 
any two consecutive polynomials Pk(x) and pk+i(x) of the sequence SRemS(p, q). In particular, 
SRemS(p, q) does not have any degree breakdown, all the e& are equal to +1, and N(q/p) is 
positive definite. 

We have, by Proposition 2.4 

q(Cj) = H n (q/p)P cn 
Since H n (q/p) is positive definite, the Cholesky algorithm gives a decomposition 

H n (q/p) = KK T 

where K 6 R nxn is lower triangular. So that we can write 

p(x) = det(zld - K- x C T p K) 

We shall remark at this point that the matrix K~ 1 C p r K is tridiagonal and symmetric. 

We get q(Cj) = KAdL where L = AdK T P CH - Then, we observe that L is a lower triangular 
matrix (since Pcw-A-d is upper triangular) and KAdL commute with Cj . Thus, we have the 
identity : 

LC^L' 1 = Ad(K- 1 CjK)Ad 

We denote by Td this tridiagonal matrix. Let (pk(x)) be the signed remainders sequence 
associated to Td as given in Proposition 1.1 (i). The first row of KAdL is proportional to 
the last row of the matrix L which is proportional to pi(x). It remains to observe that the first 
row of KAdL = q{Cp) gives exactly the coefficients of the polynomial q(x) in the canonical 
basis. Then, p\(x) = q(x). 

In summary, we have shown that, if p{x) has n simple real roots and q(x) = p'(x)/n, then 
H n (q/p) is positive definite with Cholesky decomposition H n (q/p) = KK T , and if we denote 
by q(x) the monic polynomial whose coefficients are proportional to the last row of K -1 , then 
Td(p,q) = Td(p,q). Which settle the announced duality. 

2.4 Generic case 

We turn now to the generic situation. Let p(x) and q(x) be monic polynomials of respective 
degrees n and n — 1 and such that SRemS(p, q) does not have any degree breakdown. This 
condition is equivalent to saying that all the principal minors of the Hankel matrix H n (q/p) do 
not vanish. We refer to [BPR] for this point. One way to see this is to figure out the connexion 
with the subresultants of p{x) and q{x). 

A little bit more precisely, the j-th signed subresultant coefficient of p(x) and q(x) is denoted 
by sReSj(p, q) for j = 0...n— 1. If for all j, sRes^ (p, q) ^ 0, we say that the sequence of 
subresultants is non- defective. Then, by [BPR, Corollary 8.33] and Proposition 1.1 (iv), we 
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deduce that the non-defective condition is equivalent to the fact that SRemS(p, q) has no degree 
breakdown. Moreover, from [BPR, Lemma 9.26] we know that 

(Vj € {1 . . . n}) (sRes n _>, q) = det(Hj(q/p))) . 

In conclusion, our no degree breakdown assumption means also that all the principal minors of 
the Hankel matrix H n (q/p) do not vanish. 

At this point, we may add another equivalent condition, which will be essential for the 
following. Indeed, the condition that all the principal minors of the Hankel matrix H n (q/p) 
do not vanish is also equivalent to saying that the matrix H n (q/p) admits an invertible 
LU decomposition. Namely, it exists a lower triangular matrix L with 1 entries onto the 
diagonal, and an upper invertible triangular matrix U such that H n (q/p) = LU. Moreover this 
decomposition is unique and since H n (q/p) is symmetric we may write it as H n (q/p) = LDL T 
where D is diagonal. In fact, for our purpose, we will prefer the unique decomposition 
H n (q/p) = KJK T where K is lower triangular and J is a signature matrix. 

Generalizing the previous section, we get : 

Theorem 2.5 Let p{x) and q(x) be two monic polynomials of respective degrees n and n — 1 
such that SRemS(j>, q) does not have any degree breakdown. Consider the symmetric LU- 
decomposition of the Hankel matrix H n {q/p) = KJK T , where J is a signature matrix and 
K a lower triangular matrix, and denote by q(x) the monic polynomials whose coefficients in 
the canonical basis are proportional to the last row of K^ 1 . Then, 

Td(p,q)=Td(p,q). 

Proof : We start with the companion identity : 

p(x) = det(xld - Cp) 

Next, because of Proposition 2.2(a), we notice that the matrix H n (q/p) verifies the intertwinning 
relation : 

H n (q/p)C p = C^H n (q/p) 
Then, we write the symmetric Ld-decomposition of H n (q/p) : 

H n (q/p) = KJK T 

Which gives the identity 

p(x) = det(xld - K^C^K). 

We have, by Proposition 2.4 

g(Cj) = H n {q/p)P cn = KAdL 

where 

L = AdJK T P C n 
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We observe first that L is a lower triangular matrix (since Pcw-A-d is upper triangular), and 
second that KAdL commute with Cj. Thus, we have the identity : 

LC^L- 1 = Ad(K- 1 CjK)Ad 

Proposition 1.3 gives 

Td(p, q) = Ad^^Cji^Ad 

Moreover, the first row of KAdL is proportional to the last row of the matrix L. It remains 
to observe that the first row of KAdL = q(Cj) gives exactly the coefficients of the polynomial 
q(x) in the canonical basis. Thus, by Proposition 1.3 we get 

LC^L- 1 =Td(p,q) 

Which concludes the proof. O 

Remark 2.6 Note that K _1 Cj KJ is symmetric and hence also the matrix LCj L^ 1 J , where 
J = Ad J Ad. 



3 Tridiagonal determinantal representations 
3.1 Notations 

We say that an univariate polynomial p(x) G R[x] of degree n such that p(0) / has a 
determinantal representation if 

(DR) p(x) = a det(J - Ax) 

where a £ M*, J is a signature matrix in R™ xri , and A is a symmetric matrix in R nxn (we 
obviously have a = det(J)p(0)). 

Likewise, we say that p(x) has a weak determinantal representation if 

(WDR) p(x) = a det(S - Ax) 
where a £ M*, S is symmetric invertible and A is symmetric. 

Of course the existence of (DR) is obvious for univariate polynomials, but we will focus on 
the problem of effectivity. Namely, we want an algorithm (say of polynomial complexity with 
respect to the coefficients and the degree of p{x)) which produces the representation. Typically, 
we do want to avoid the use of the roots of p(x). 

One result in that direction can be found in [Qz2] (which is inspired from [Fi]). It uses arrow 
matrices as a "model", whereas in the present article we make use of tridiagonal matrices. 

When all the roots of p(x) are real, the effective existence of determinantal representation 
for univariate real polynomials exists even if we add the condition that J = Id. It has 
been discussed in several places, although not exactly with the determinantal representation 
formulation. Indeed, in place of looking for DR we may consider the equivalent problem of the 
research of a symmetric matrix whose characteristic polynomial is given. Indeed, if the size of 
the matrix A is equal to the degree n of the polynomial, the condition 

p(x) = det(Id — xA) 
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is equivalent to 

p*(x) = det(xld - A) 

where p*(x) is the reciprocal polynomial of p(x). In [Fi], arrow matrices are used to answer 
this last problem. On the other hand, the Routh-Lanczos algorithm (which can be viewed as 
Proposition 1.1) gives also an answer, using tridiagonal model. Note that the problem may also 
be reformulated as a structured Jacobi inverse problem (confer [EP] for a survey). 

In the following, we generalize the tridiagonal model to any polynomial p(x), possibly having 
non real roots. Doing that, general signature matrices J appear, whose entries depend on the 
number of real roots of p{x). 

3.2 Over a general field 

A lot of identities in Section 2 are still valid over a general field k. For instance, if p(x) and 
q(x) are monic polynomials of respective degrees n — 1 and n, we may still associate the Hankel 
matrix H(q/p) = (sj + j)o<jj< n -i G k nxn defined by the identity 

P ( X ) Zj, x j+i 

Then, we have the following : 

Theorem 3.1 Let p{x) 6 k[x], q(x) 6 k[x] be two monic polynomials of respective degrees n 
and n — 1, and set H = H(q/p). Then, the matrix C^H is symmetric and we have the WDR : 

det(H) x p(x) = det(xH - C^H) 

Moreover, if H admits the LU- decomposition H = KDK T where K G k nxn is lower triangular 
with entries 1 onto the diagonal and D € k nxn a diagonal matrix, then we have : 

p(x) = det(xD - Td) 

where Td = K~ 1 CpKD is a tridiagonal symmetric matrix. 

Proof : We exactly follow the proof of Theorem 2.5. O 

Note that the condition for H to be invertible is equivalent to the fact that the polynomials 
p(x) and q(x) are coprime, since we have 

rk(Bez(g,p)) = deg (p) - deg (gcd(p, q)). 

To see this, we may refer to the first assertion of [BPR, Theorem 9.4] whose proof is valid over 
any field. 

The WDR of Theorem 3.1 has the advantage that the considered matrices have entries in 
the ring generated by the coefficients of the polynomial p(x). This point is not satisfied in the 
methods proposed in [Qz2] or in the Routh-Lanczos algorithm. 

In fact, the use of Hankel matrices satisfying the intertwinning relation seems to be more 
convenient since we are able to "stop the algorithm at an earlier stage", namely before having 
to compute a square root of the matrix H (or of the matrix D). 
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Of course, at the time we want to derive a DR, then we have to add some conditions on the 
field k, for instance we shall work over an ordered field where square roots of positive elements 
exist. 

To end the section, we may summarize that, for a given polynomial p(x), we have an obvious 
but non effective (i.e. using factorization) DR with entries in the splitting field of p(x) over k, to 
compare with an effective WDR given by Theorem 3.1 where entries are in the field generated 
by the coefficients of p(x). 

3.3 Symmetric tridiagonal representation and real roots counting 

If p{x) and r(x) are two real polynomials, we recall the number known as the Tarski Query : 

TaQ(r,p) = #{x G M | p(x) = A r(x) > 0} - #{x G M | p(x) = A r(x) < 0}. 

We also recall the definition of the Permanences minus variations number of a given sequence 
of signs v= (i/i, . . . ,i/ fc ) : 

fc-i 

PmV(^) = ^2 vm+i- 

i=l 

We summarize, from [BPR, Theorem 4.32, Proposition 9.25, Corollary 9.8] some useful 
properties of these numbers, 

Proposition 3.2 Let p(x) and q(x) be two monic polynomials of respective degrees n andn—1, 
and such that the sequence SRemS(p, q) has no degree breakdown. Let r(x) be another polynomial 
such that q(x) is the remainder of p'(x)r(x) modulo p(x). Then, 

PmV(i/) = sgn(Bez(p,g)) = sgn(H n (q/p)) = TaQ(r,p) 

where v is the sequence of signs of the leading coefficients in the signed remainders sequence 
SRemSO, q). 

We come now to our main result about real roots counting : 

Theorem 3.3 Let Td G R nxn be a tridiagonal symmetric matrix with non-singular first 
principal diagonals. Let also p(x) G R[x] be a real polynomial with no multiple root and such 
that 

p(x) = det(J) det(xJ - Td), 

where J is a signature matrix whose last entry onto the diagonal is +1. 
Then, the number of real roots of p(x) is greater than sgn(J). 

Proof : We have 

p(x) = det(xld n - Td x J) 

and we set 

q(x) = det (zld n _i - (Td x J) n-1 ) . 

The matrix Td x J is still tridiagonal with non-singular first principal diagonals. Then, we 
consider the sequence SRemS(p, q) and denote by v the associated sequence of signs of leading 
coefficients. 
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Since gcd(p,p') = 1, we set r(x) to be the unique polynomial of degree < n such that 



r = — mod p. 
pi 



Then, 



p'r = q mod p 

and from Proposition 3.2, we get : 

PmV(^) = TaQ(r,p) < #{x £ R | p(x) = 0} 

Let us introduce some notations at this step. First, let Td = Td(a,/3, 7), next denote by 
e(a) the sign in { — 1, +1} of the non zero real number a, and finally let 



J 



/ 0n- 



\ 



01 



\ 



1 / 



Then, we can write 



where 



p(x) = det (xld n - P(Td x J)P~ l ) 
( {0 n -i...0{) x (e( 7n _i) . . . e( 7 i)) 



P = 



01 x e(7i) 



V 



1 / 



We note in fact that P(Td x J)P _1 = Td(p, g). Indeed, all the coefficients onto the first 
lower principal diagonal are positive. Moreover, all the coefficients onto the first upper principal 
diagonal are given by the sequence 

(9 n -i x 6 n -2, ■ ■ ■ , O2 x 6\,9i). 

We deduce from Proposition 1.1 (iv) that the sequence of signs of leading coefficients in the 
signed remainders sequence SRemS(p, q) is the following : 



Thus 



and we are done. 



V = {On-! X . . . X 0i, . . . , 6 2 X #1, 01, 1, 1). 
n-1 

PmVH = l + ^0 fc = S gn(J) 



fc=i 



□ 



Another way, maybe less constructive, to prove the result is to use the duality of Theorem 
2.5. Indeed, replacing as in the previous proof the matrix Td x J with P(Td x J)P _1 , we write 
the identity 



T r -1 



Td x J = LCp L 
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Then, by duality, we have 

LC^L- 1 = AdK^CjKJ'Ad 
where we have set the LU-decomposition 

H n (q/p) = KJ'K T . 

Let us introduce J' = AdJ'Ad ; we get : 

(LCjL _1 j) x (J J') = AdK~ 1 C p KJ'Ad 

We remark that the matrices LC p r L~ 1 J and K~ l C p KJ' are both tridiagonal and symmetric 
with non-singular principal diagonals, so that we necessarily have 

J J' = ±Id. 

Notice that by assumption the last coefficient of J is +1 and that the first coefficient of J' is 
always +1 (since it is the leading coefficient of Thus 

JJi = Id. 

Then, we may conclude by Proposition 3.2. 

An alternative way to make use of this computation is to say that we get another proof of 
the equality 

PmV(^) = sgn(Bez(p, q)) 
which appears in the sequence of identities 

sgn(Bez(p, q)) = sgn(H n (q/p)) = sgn(j') = sgn(J) = PmV(z^) = TaQ(r,_p). 

Remark 3.4 It is possible to extend Theorem 3.3 in the case where principal diagonals of 
Td = Td(a, (3, (3) are singular. Namely, for all k such that [3k = 0, we have to assume that 
the corresponding k-th entry onto the diagonal of J is equal to +1. Then, we get that the 
number of real roots of p(x), counted with multiplicity, is greater than sgn(J). 

To see this, it suffices to note that the polynomial defined by p{x) = det( J) det(x J — Td) 
factorizes through 

p(x) = det(Ji) det(xJi - Td fc ) x det(J 2 ) det(xJ 2 - Td„_ fc ) 

Moreover, the matrices Tdk and Td n _fc remain tridiagonal symmetric and J\, J 2 remain 
signature matrices. If we denote by the usual direct sum of matrices, we have J = J\ J 2 
and Td = Tdjfc0Td n _ fe . 

Thus, we may proceed by induction on the degree of p{x). 

Before stating the converse property of Theorem 3.3, we establish a genericity lemma. 

Lemma 3.5 Let p(x) be a monic polynomial of degree n with only single roots and q(x) = 
x n-i _|_ i lX n-i _|_____|_ b n _ 1- Then, the set of all (n — l)-tuples (b\, . . . , 6 n _i) 6 ]R ra_1 such that 
there is an integer k G {1, . . . , n} satisfying det(Hk(q/p)) = 0, is a proper subvariety o/M n_1 . 
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Proof: We only have to show that for all k, det(Hk(q/p)), viewed as a polynomial in the 
variables b\, ... , 6 n _i, is not the zero polynomial. 
Let H n (q/p) = (s i+j )o<i,j< n -i where 



p(x) 



3=0 




and denote by a>i, . . . , a n the set of all (possibly complex) roots of p(x). Then, 



n 



P'{aci) 



Let us introduce the real numbers defined as 



n 




u i = Yl 



P'{aci) 



i=l 



We obviously have Uj = whenever j < n — 2 and also « n -i 
that we deduce : 



= 1 (look at lim. 



x j q(x) 



). So 



si = h + u n 
and more generally 
k (Vj G {1, ... ,2n — 2}) (sj = 6j + bj-\u n + . . . + biu n+j - 2 + u n +j-i) 

Then, it becomes clear that Hk + i(q/p) ^ for any k such that k < L^y^J = r > since 
S2k G R[6i, . . . , 62/c] has degree 1 in the variable &2fc and so is the case for Hk+i(q/p). 

Next, for r < k < n, we develop the determinant Hk{q/p) successively according to the first 
columns, and we remark that its degree in the variable 6 n _i is equal to 2k — n (with leading 
coefficient equal to —1). This concludes the proof. O 

In other words, the Lemma says that the condition 



is generic with respect to the space of coefficients of the polynomial q(x). Because of the relations 
between coefficients and roots, the condition is also generic with respect to the (possibly complex) 
roots of the polynomial q(x). 

Here is our converse statement about real roots counting : 

Theorem 3.6 Letp{x) be a monic polynomial of degree n which has exactly s real roots counted 
with multiplicity. We can find effectively a generic family of symmetric tridiagonal matrices Td 
and signature matrices J with sgn( J) = s, and such that 



Proof: If p{x) has multiple roots, then we may factorize it by gcd(p,j/) and use the 
multiplicative property of the determinant to argue by induction on the degree. Now, we assume 
that p{x) has only simple roots. 



(VfcG {!,..., n}) (det(H k (q/p)) = 0) 



p(x) = det(J) x det(xJ-Td). 
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We take for q{x) any monic polynomials of degree n — 1 which has exactly s — 1 real roots 
interlacing those of p{x). Namely, if we denote by x\ < . . . < x s all the real roots of p(x) and 
by yi < ... < Us-i ah the real roots of q(x), we ask that x\ < y\ < x\ < yi < . . . < y s -i < x s . 

Let r{x) be the unique polynomial of degree < n such that r(x) = mod p{x) (since p'(x) 
is invertible modulo p(x)). 

From p'r = q mod p and p'(xi) = q(xi) for all real root Xj of p(x), we get 

TaQ(r,p) = s = #{x G E | p(x) = 0} 

At this point, we need that q(x) satisfies another hypothesis : that is SRemS(p, q) shall 
not have any degree breakdown, or equivalently that H(q/p) shall admit a L?7-decomposition 
H n (q/p) = KJK T . According to Lemma 3.5, this hypothesis is generically satisfied, although 
it may not be always satisfied for the natural candidate q(x) =p'(x)/n. 

Then, we get from Theorem 2.5 

p{x) = det(xJ - K^C^KJ) 

where Td = K~ 1 CpKJ is tridiagonal symmetric and J is a signature matrix. 

By the proof of Proposition 3.3, we get moreover that 

sgn(J) = TaQ(r,p) = S gn(H n (q/p)). 

This concludes the proof since TaQ(r, p) = s. O 

Remark 3.7 (i) The choice of such polynomials q(x) with the interlacing roots property need 
to count and localize the real roots ofp(x) . It can be done via Sturm sequences for instance. 

(ii) Although the polynomial q(x) = p'(x)/n has not necessarily the interlacing property in 
general, it is the case when all the roots of p(x) are real and simple. Moreover, in this 
case, the interlacing roots condition is equivalent to the no degree breakdown condition. 
Indeed, TaQ(p'g mod p,p) = n if and only if p'(x) and q(x) have same signs at each root 
ofp(x). 

4 Some worked examples 

In order to get lighter formulas in our examples, we decide to get rid off denominators. That 
is why we replace signature matrices by only non-singular diagonal matrices. If one wants to 
deduce formulas with signature matrices, it suffices to normalize. 

1) Let p(x) = x 3 + sx + t with s / 0, and q(x) = p'(x) = 3x 2 + s. Let us introduce the 
discriminant of p(x) as A = — 4s 3 — 27t 2 . Consider the decomposition of the Hankel matrix 

/ 3 -2s \ 
H(q/p) = -2s -3t = KJK T 
\ -2s -3t 2s 2 / 
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where 



and 




We recover the well-known fact that p(x) has three distinct real roots if and only if s < 
and A > 0, which obviously reduces to the single condition A > 0. 

Then, we have the determinantal representation 



A x p(x) = det(xJ - Td) 



where 




2) Consider the polynomial p(x) = x 5 — 5x 3 + 4x, which in fact factorizes through p(x) 
x(x - l)(x + l)(x - 2)(x + 2). Let q(x) = p'(x)/5. We have 



N(q/p) 



( 5 





10 





34 \ 





10 





34 





10 





34 





130 





34 





130 





V 34 





130 





514 / 



Td 



/ V2 

V2 





y 








o \ 




/3b o 



p{x) = det(xId-Td). 

In order to get some parametrized identities, let us introduce the following family of 
polynomials 

3\ / 1\ / 1' 



q a {x) = {x - a) ( x + - ] \ x +2) \ x ~ 2 
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We write the LU-decomposition 



H(q a /p) 



\ 



1 


3 
2 


- a 


—3a 
2 


+ 19 
^ 4 


f-a 


—3a 
2 


^ 4 


57 
8 


19a 
4 


-3a , 19 
2 4 


57 
8 


19a 

4 


-57a 
8 


+ 79 
+ 4 


57 19a 


—57a 


+ 79 


237 


79a 


8 4 


8 


^ 4 


8 


4 


-57a , 79 
8 ~r 4 


237 
8 


79a 
4 


-237a 
8 


, 319 
~r 4 



57 
8 


19a 
4 


-57a 
8 


+ 79 
^ 4 


\ 


-57a 


f ^ 


237 


79a 




8 


r 4 


8 


4 




237 
8 


79a 
4 


-237a 
8 


, 319 
"T" 4 




-237a 


, 319 


957 


319a 




8 


r 4 


8 


4 




957 

8 


319a 
4 


-957a 
8 


i 1279 
4 


/ 



K a J a K T a 



where the associated "signature" matrix J a is equal to 
/ 1 



■i(o+l)(2o-5) 



/15\ (2a-l)(4a 2 -q-15) 
I 16/ (a + l)(2a-5) 



/ 45 \ 48a 4 - 16a 3 -216a 2 +58a+ IPS 
I 128 / 



V 



(2a-l)(4a 2 -a-15) 



/315\ (a + 2)(a + l)a(a-l)(a-2) 
IS/ 48a 4 -16a a -216a 2 +58a + 105 



The condition for H a (q/p) to be positive definite is equivalent to having only positive 
coefficients onto the diagonal of J a . 

First, it yields J a (2,2) > 0, which means that a g] — 1, §[. Then, we add the condition 
J a (3,3) > which means that a g]|,2, 06..[. Then, we add the condition J a (4, 4) > 
which means that a e]0, 9.., 2, 00..[. And finally, we add the condition J a (5, 5) > 0, which 
means that a 2[ and gives exactly the interlacing property for the polynomial q a (x). 



For instance, with a = | we get p(x) = det (xld — Td3^ where : 



/ 



Tda 

2 















/35 q 



\ 
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